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1. Introduction 

Chiral Perturbation Theory (ChPT) [|], ^ describes the low energy interactions of pions, which 
are the pseudo-Goldstone bosons associated to the spontaneous chiral symmetry breaking of QCD, 
and hence, the relevant degrees of freedom at low energies. ChPT is built as the most general low 
energy expansion in terms of the pion momenta and mass that is compatible with the QCD symme- 
tries. The details of the underlying dynamics at higher energies are encoded in a set of parameters, 
known as low energy constants (LECs), that appear at the different orders in the expansion, and, 
once renormalized, absorb the loop divergences present at that order. Since perturbative QCD can- 
not be applied at very low energies, it is particularly difficult to obtain the values of these LECs 
from first principles and, with few exceptions, the LECs have been determined best from the com- 
parison with experiment [|[ ^ ^, ^. There are determinations of the LECs from lattice QCD (we 
refer to for a recent compilation), and it is also possible to obtain positivity constraints for the 
LECs coming from axiomatic field theory 

In the ChPT nn scattering amplitude only certain combinations of the LECs appear to a given 
order of the calculation. To leading order, 0{p^), there are no LECs and only the pion mass decay 
constant M^i and appear. To next to leading order (NLO), or 0{p'^), which corresponds to a one- 
loop calculation, only four LECs, called appear in the amplitude. At next to next to leading 
order (NNLO), or 0{p^), which corresponds to a two-loop calculation, only six independent terms 
appear [||], multiplied by the corresponding combination of LECs, denoted bi, i = 1...6. 

Here we report our work [||] on the determination of the LECs that appear in the KK scattering 
amplitude. We obtain these LECs from fits to the coefficients of the momentum expansion of the 
amplitude around threshold, known as threshold parameters, which can be calculated within ChPT 
and expressed in terms of the LECs. The experimental value of the threshold parameters we use is 
the one obtained from a dispersive data analysis in [TH], which includes the very precise and reliable 



results on K14 decays from the NA48/2 collaboration [12|. We also provide here the calculation of 



the third order coefficients of the threshold expansion, which was not performed in [p^. 

2. Threshold parameters 

The amplitude for tttt scattering is customarily decomposed in terms of partial waves t\, of 
definite isospin / and angular momentum t. t\{s) = /Ij {s,t ,u)Pe{cos 6) d{cos 6), 6 being 
the scattering angle. Pi the Legendre polynomials, s,t,u the usual Mandelstam variables and T the 
amplitude. With this normalization, the threshold expansion can be written as: 

■^Re4{s)=p^^(afj + beip^ + ^ceip^ + ..}j, (2.1) 

where the aa are usually called scattering lengths, the ba slope parameters, the cei shape parame- 
ters, and all of them, generically, threshold parameters. 

The use of sum rules to obtain the values of threshold parameters is a well established tech- 
nique (pi [nil that we will also use here. Our experimental determination of the threshold parameters 



is performed using the parametrizations of [10|, which were obtained by highly constraining data 
fits to satisfy three sets of dispersion relations within uncertainties. In [p^], the values of the a 
and b parameters up to F waves were provided. With the aim of minimizing the uncertainties, they 
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were obtained from sum rules, with the only exception of the Saso + 2^52 combination, which is 
orthogonal to the one appearing in the Olsson sum rule (note the spectroscopic notation, where the 
£ = 0, 1,2,3... are denoted S,P,D,F...). 

Here, we also provide the calculation of the third order coefficients c of the threshold expan- 
sion, which adds five more observables for the fit. For the c parameters with £ > 0, we use the 
Froissart-Gribov sum rules: 
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where F'{s,t) = T'{s,t)/A7Z^ and the primes denote the derivative with respect to cos0. This 
formula allow us to calculate the c parameters for the P, DO, D2 and F waves. For the S waves we 
provide two new sum rules, and also one for cp, in order to reduce its error: 
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The derivation is similar to that of the sum rules for bp, bso and bsi obtained in []13[]. They corre- 
spond to the threshold limit, taken from above, of the second derivative of a forward dispersion rela- 
tion for the F'^^, F^'^ and F^ amplitudes, respectively. Let us recall that F^^ = F'^'^ /2 + F'^^ /2 
whereas F^ = 2F'^^ /3 + F'^^/3. For a list of the resulting values, we refer the reader to Table II 
of our original work 



3. 0(/)fits 

We start by fitting the 0{p^) LECs by using the one-loop expression of the threshold parame- 
ters, which will help us check the stability of the LECs values and the need for higher orders. We 
actually fit the It parameters, which are basically the /[(jU) at the n = scale and normalized so 
that they have values of order one [Q]. Note, however, that It, and l/[ only appear through the quark 
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Fit to 


h 


h 


h 


k 


X^d.o.f. 


as,bs,ap 


l.lil.O 


5.1±0.7 


-1±8 


7.1±0.7 


0.23 


ao 


-1.75±0.22 


5.91±0.10 









cs 


-2.4±0.9 


4.8±0.4 









as,bs,ap,aD,cs,bp 


-2.06±0.14 


5.97±0.07 


-5±8 


7.1±0.6 


7.9 


as,bs,ap,aD,cs,bp, using /o 


-1.06±0.11 


4.6±0.9 


0±6 


5.0±0.3 


7.06 


Estimate 0{p'^) 


-1.5±0.5 


5.3±0.7 


-3±7 


6.0±1.2 





Table 1: 0{p'^) fits to different sets of threshold parameters containing polynomial 0[p'^) contributions. We 
observe that a precise description of the observables is not possible at one loop. Anyway, we provide an 
estimate of how much one should enlarge the uncertainties of the LECs if, for simplicity, one still insists in 
using the one-loop formalism. 



mass dependence of Mji and fji, respectively, and therefore we cannot expect much sensitivity to 
these two parameters from fits to the coefficients of the momentum expansion of amplitudes. In 
addition, since the LECs only appear in the polynomial part of the partial waves, which at one loop 
is of 0{p^), only ten observables carry any dependence on the LECs: a^o), asi)^ ^so> ^52> bp, 
cso> cs2, ^DQ and a^. The rest of the coefficients multiply powers of the momentum higher than 
and thus, do not receive a contribution from the 0{p^) LECs. 

In Table I we show the results of our fits. First, we have fitted only the observables whose 
leading contribution is of 0{p^), since these might be more stable under the higher order correc- 
tions. The fit comes out with relatively low jd.o.f .. Next we present two determinations of 
l\ and I2, which can be fixed using only either am and am, or C50 and 0^2- It is evident that the 
resulting values from those fits are incompatible. The incompatibility is even worse when fitting 
simultaneously the ten observables that depend on Ij to 0{p'^), where we obtain a high x^/d.o.f. 
value. Finally, the effect of higher order corrections has been studied by fitting to the one-loop 
amplitude but replacing fj^ by /o in the 0{p'^) terms, since the two expressions only differ in higher 
order contributions. This we show in row 5 of Table || The x^/d.o.f is somewhat lower, but the 
values of the LECs come out rather different from the previous calculation. 

These results imply that, to the present level of precision, the one-loop ChPT formalism is not 
enough and calls for higher order corrections. If one still wants to use this simpler version instead 
of the full two-loop amplitude one can include the effect of higher orders into a systematic uncer- 
tainty of the LECs. We propose to take the weighted average of the two previous fits, including 
a systematic uncertainty to cover the LECs values of both fits. In Table IV of the original refer- 
ence [^], we compare the resulting threshold parameters obtained using this averaged set with the 
experimental values. Thanks to the larger uncertainty, the threshold parameters obtained are com- 
patible within errors with the experimental values, except for ^50 and bp, which differ by more than 
three and two standard deviations respectively. It is worth noting also that the positivity constraints 
obtained from first principles are perfectly satisfied by this set of LECs, even in the worst case 
scenario. 
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4. 0{p^) fits 

As commented in the introduction, the two-loop nn scattering amplitude can be recast in terms 
of six independent terms multiplied by their corresponding low energy constants hi. In turn, these bj 
can be rewritten in terms of the four 0{p^) LECs and six combinations r, of 0{p^) LECs [14]. The 
difference in the amplitude using one way or the other is 0{p^). However, despite increasing the 
number of parameters to ten, the 0{p^) amplitude still provides just six independent structures. As 
a consequence, the fits in terms of It and r, are much more unstable, and can even lead to spurious 
solutions. For this reason we focus on the fits in terms of bi, and refer to the appendix for a 
study of the It, r,- fits. 

We first fit the ten threshold parameters used in the previous section because, having a non-zero 
0{p^) polynomial contribution, we expect these to be more stable under higher order corrections. 
In the first row of Table ||we show the resulting bj, which describes fairly well the fitted observables 
with a /d.o.f. = 1.2. However, when fitting all 18 observables, we obtain somewhat different 
LECs (see the second row of Table ^) and the jd.o.f . comes out rather poor. We have noticed 
that cp alone contributes almost to one third of the total 'X^ . This might indicate that cp receives 
important higher order contributions that are not being taken into account in the 0{p^) calculation. 
Once again we obtain a crude estimate of the size of higher order ChPT corrections, by changing 
fn by /o in the last term of the expansion, cp suffers indeed the largest change, by almost 80%. 

Thus, we proceed to fit again all threshold parameters except cp. The result is shown in the 
third row of Table ^ The fit quality improves sizably, but we still get a high /d.o.f. = 2.9, 
which indicates that the two-loop calculation may not be enough to describe even the remaining 
threshold parameters with their current level of precision. Again, we see the effect of higher order 
corrections by making a fit replacing by /o in the 0{p^) terms. We show the results in the 
fourth row of Table ^. Surprisingly, we now obtain a good x^ M-o.f. = 1.0 and all LECs are less 
than two standard deviations away from those obtained by fitting only the threshold parameters 
with an 0{p'^) polynomial part. We conclude that, by excluding cp, the two-loop fit can give an 
acceptable description of the rest of threshold parameters. For this reason, we have once more 
made a weighted average of the two fits (the one using and the one using /o) adding systematic 
uncertainties to cover both sets. This we show in the fifth row of Table ^ where we can see that 
they are also quite compatible with previous determinations in the literature [||]. Also, this LECs 
satisfy again the axiomatic constraints 



5. Summary 

We have reported a work [||] on the determination the low energy constants of SU (2) Chiral 
Perturbation Theory (ChPT) at one and two loops by fitting to the threshold parameters obtained 



from sum rules using a recent and precise dispersive analysis of data [10], together with six addi- 
tional observables that we have studied here. 

We have checked that the one-loop formalism is clearly insufficient to accommodate the 
present level of precision. The x} l^-o.f . improves remarkably when using the two-loop expan- 
sion, although it is still not sufficient to get a good quality fit. This suggests that even higher order 
ChPT contributions may still be required to describe all these observables simultaneously. 
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Fit to 


bi 


h 


h 


b4 


^5 


b6 


X 
d.o.f. 


as,bs,ap,aD,cs,bp 


-14±4 


14.6±1.2 


-0.29±0.05 


0.76±0.02 


O.lil.l 


2.2±0.2 


1.2 


AU 


-2±3 


14.2±1.0 


-0.39±0.04 


0.746±0.013 


3.1±0.3 


2.58±0.12 


5.2 


All but cp 


-6±3 


15.9±1.0 


-0.36±0.04 


0.753±0.013 


2.2±0.4 


2.44±0.12 


2.9 


All but Cp, using /o 


-12±3 


13.9±0.9 


-0.30±0.04 


0.726±0.013 


1.0±0.3 


1.93±0.08 


1.04 


Estimate 0{p^) 


-10.5±5.1 


14.5±1.8 


-0.31±0.06 


0.73±0.02 


1.3±1.0 


2.1±0.4 




Ref. [| 


-12.4±1.6 


11.8±0.6 


-0.33±0.07 


0.74±0.01 


3.6±0.4 


2.35±0.02 





Table 2: 0{p'') fits. In the first row we only fit to observables containing polynomial 0{p'^) contributions. 
Next we show the fit to all the threshold parameters obtained in this work. The quality is rather poor, but 
most of the disagreement is caused by cp. When this observable is omitted, the resulting fits are much better, 
specially when using /q instead of fji in the last term of the ChPT expansion. We provide an estimate of the 
LECs uncertainties from the fits to all observables except cp, as a weighted average of the fits using /o or fj^. 
The resulting hj parameters are very consistent with previous determinations, listed in the last row. 
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